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Abstract 

There is vast empirical evidence that given a set of assumptions on 
the real- world dynamics of an asset, the European options on this asset 
are not efficiently priced in options markets, giving rise to arbitrage 
opportunities. We study these opportunities in a generic stochastic 
volatility model and exhibit the strategies which maximize the arbi- 
trage profit. In the case when the misspecified dynamics is a classical 
Black-Scholes one, we give a new interpretation of the butterfly and 
risk reversal contracts in terms of their performance for volatility ar- 
bitrage. Our results are illustrated by a numerical example including 
transaction costs. 
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1 Introduction 



It has been observed by several authors [TJ [3J [12] that given a set of assump- 
tions on the real-world dynamics of the underlying, the European options 
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Figure 1: Historical evolution of the VIX index (implied volatility of options 
on the S&P 500 index with 1 to 2 month to expiry, averaged over strikes, see 
[10J ) compared to the historical volatility of the S&P 500 index. The fact 
that the implied volatility is consistently above the historical one by several 
percentage points suggests a possibility of mispricing. 

on this underlying are not efficiently priced in options markets. Important 
discrepancies between the implied volatility and historical volatility levels, as 
illustrated in Figure [Tj as well as substantial differences between historical 
and option-based measures of skewness and kurtosis [3] have been docu- 
mented. These discrepancies could be explained by systematic mispricings / 
model misspecification in option markets, leading to potential arbitrage op- 
portunities. The aim of this paper is to quantify these opportunities within a 
generic stochastic volatility framework, and to construct the strategies max- 
imizing the gain. The arbitrage opportunities analyzed in this paper can be 
called statistical arbitrage opportunities, because their presence / absence 
depends on the statistical model for the dynamics of the underlying asset. 
Contrary to model independent arbitrages, such as violation of the call-put 
parity, a statistical arbitrage only exists in relation to the particular pricing 
model. 

The issue of quantifying the gain/loss from trading with a misspecified 

1 There exist many alternative explanations for why the implied volatilities are con- 
sistently higher than historical volatilities, such as, price discontinuity [3], market crash 
fears [5] and liquidity effects such as transaction costs [22l [T9l [8] [9] , inability to trade in 
continuous time [6] and market microstructure effects (24]. The literature is too vast to 
cite even the principal contributions here 
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model has so far mainly been studied in the case of the Black-Scholes model 
with misspecified volatility [13, 25J. In this paper we go one step further, and 
analyze the effects of misspecification of the volatility itself, the volatility of 
volatility and of the correlation between the underlying asset and the volatil- 
ity in a stochastic volatility model. Since these parameters may be observed 
from a single trajectory of the underlying in an almost sure way, their mis- 
specification leads, in principle, to an arbitrage opportunity. The questions 
are whether this opportunity can be realized with a feasible strategy, and 
how to construct a strategy maximizing the arbitrage gain under suitable 
conditions guaranteeing the well-posedness of the optimization problem. 

While the issue of consistency between real-world and risk-neutral prob- 
ability measures has been given a rigorous treatment in several papers [TJ [3j 
H2] , the corresponding arbitrage trading strategies are usually constructed in 
an ad- hoc manner [TJ [T71 [18]. For instance, when the risk- neutral skewness 
is greater than the historical one (which roughly corresponds to correlation 
misspecification in a stochastic volatility model), [lj suggest a strategy con- 
sisting in buying all OTM puts and selling all OTM calls. Similarly, if the 
risk-neutral kurtosis is greater than the historical one, the trading strategy 
consists in selling far OTM and ATM options while simultaneously buying 
near OTM options. In this paper we determine exactly which options must 
be bought and sold to maximize arbitrage gains, depending on model param- 
eters. 

Our second objective is to analyze commonly used option trading strate- 
gies, such as butterflies and risk reversals, and provide a new interpretation 
of these structures in terms of their performance for volatility arbitrage. A 
butterfly (BF) is a common strategy in FX option trading, which consists 
in buying an out of the money call and an out of the money put with the 
same delta value (in absolute terms) and selling a certain number of at the 
money calls/puts. A risk reversal (RR) is a strategy consisting in buying an 
out of the money call and selling an out of the money put with the same 
delta value (in absolute terms). The financial engineering folklore goes that 
"butterflies can be used to arbitrage misspecified volatility of volatility" and 
"risk reversals can be used to arbitrage misspecified correlation". In section HI 
we study these strategies and discuss their optimality for volatility trading. 

During the last decade we have witnessed the appearence of a large spec- 
trum of new products specifically designed for volatility and correlation trad- 
ing, such as variance, volatility and correlation swaps. However, in most mar- 
kets, European options continue to be much more liquid than exotic volatility 
products and still constitute the most commonly used tool for volatility arbi- 
trage. In this paper we therefore concentrate on arbitrage strategies involving 
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only the underlying asset and liquid European options. 

The rest of the paper is structured as follows. In Section [21 we introduce 
the generic misspecified stochastic volatility model. Section [3] defines the 
admissible trading strategies and establishes the form of the optimal arbitrage 
portfolio. Section H] is dedicated to the special case when the misspecified 
model is the constant volatility Black-Scholes model. This allows us to give 
a new interpretation of butterfly spreads and risk reversals in terms of their 
suitability for volatility arbitrage. Section [5] presents a simulation study of 
the performance of the optimal arbitrage strategies in the framework of the 
SABR stochastic volatility model [T5]. The purpose of this last section is not 
to prove the efficiency of our strategies in real markets but simply to provide 
an illustration using simulated data. A comprehensive empirical study using 
real market prices is left for further research. 



2 A misspecified stochastic volatility framework 

We start with a filtered probabiblity space (O, J 7 , P, (J r t)t>o) and consider a 
financial market where there is a risky asset S, a risk-free asset and a certain 
number of European options on S. We assume that the interest rate is zero, 
and that the price of the risky asset S satisfies the stochastic differential 
equation 

dSt/St = tHdt + a t y/l - pJdWl + a t p t dW 2 (1) 

where p, a and p G [—1,1] are adapted processes such that 
t 

(1 + S 2 S ){\ + p 2 s + cr 2 s )ds < oo a.s. for all t, 

and (W /1 ,W /2 ) is a standard 2-dimensional Brownian motion. This integra- 
bility condition implies in particular that the stock price process never hits 
zero P-a.s. 

To account for a possible misspecification of the instantaneous volatility, 
we introduce the process a t , which represents the instantaneous volatility 
used by the option's market for all pricing purposes. In particular, it is 
the implied volatility of very short-term at the money options, and in the 
sequel we call a the instantaneous implied volatility process. We assume that 
<r t = &(Yt), where Y is a stochastic process with dynamics 

dY t = a t dt + b t dW t 2 , (2) 
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where at and bt > are adapted processes such that 

t 

(al + b 2 s )ds < oo a.s. for all t, 

and a : 1R — > (0, oo) is a continuously differentiable Lipschitz function with 
< a < a(y) < o < oo and a'(y) > for all y G R; 

Further, to account for possible misspecification of the volatility of volatil- 
ity b and of the correlation p, we assume that there exists another probability 
measure Q, called market or pricing probability, not necessarily equivalent 
to P, such that all options on S are priced in the market as if they were 
martingales under Q. The measure Q corresponds to the pricing rule used 
by the market participants, which may be inconsistent with the real-world 
dynamics of the underlying asset (meaning that Q is not necessarily abso- 
lutely continuous with respect to P). Under Q, the underlying asset and its 
volatility form a 2-dimensional Markovian diffusion: 

dS t /S t = ~a{Y t )^/l-~p\Y u t)dW} + a(Y t )p(Y t ,t)dW? (3) 
dY t = a(Y t , t)dt + b(Y t , t)dW?, (4) 

where the coefficients a, b and p are deterministic functions and (W 1 ,^ 2 ) 
is a standard 2-dimensional Brownian motion under Q. Since a is bounded, 
the stock price process never hits zero Q-a.s. 

The following assumptions on the coefficients of (J3j) — (SJ) will be used 
throughout the paper: 

i) There exists e > such that min(l — p(y, t) 2 , b(y, t)) > e for all (y, t) G 
R x [0,T]. 

ii) The functions d(y, t), b(y, t), p(y, t) are twice differentiable with respect 
to y; these coefficients as well as their first and second derivatives with 
respect to y are bounded and Holder continuous in y, t. 

iii) The function a is twice differentiable; this function as well as its first 
and second derivative is bounded and Holder continuous. 

We suppose that a continuum of European options (indifferently calls or 
puts) for all strikes and at least one maturity is quoted in the market. The 
price of an option with maturity date T and pay-off H(St) can be expressed 
as a deterministic function of St, Y t and t: 

P(S t ,Y t ,t) = E Q [H(S T )\^t]. 
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Using standard methods (see e.g. [15]) one can show that under our as- 
sumptions, for every such option, the pricing function P belongs to the class 
C 2 ' 2 ' 1 ((0, oo) x R x [0,T)) and satisfies the PDE 

dP 

~a— + CP = 0, (5) 

where we define 

~ df S 2 a(y) 2 d 2 f , b 2 d 2 f ~ d 2 f 



dt 2 dS 2 2 dy 2 vy/ r dSdy' 
In addition (see [26]), the price of any such European option satisfies 
OP 

— >0, V(S, y, t) G (0, oo) xlx [0, T). (6) 
dy 

We shall use the following decay property of the derivatives of call and 
put prices (see Appendix A for the proof). 

Lemma 2.1. Let P be the price of a call or a put option with strike K and 
maturity date T . Then 

Um aP(Sy,t) = ^aP(Sy,t) 
k^+oo dy k^o dy 

d 2 P(S,y,t) d 2 P(S iy ,t) 
hm — = hm — = 0, 

k-^+oo dy k-*o dy 

d 2 P(S,y,t) d 2 P(S,y,t) 
hm — = hm — = 0, 

K-^+oo dS 2 k^o dS 2 

d 2 P(S,y,t) d 2 P(S } y,t) 

and lim = lim = 

k^+oo dSdy k^o dSdy 

for all (y,t) Glx [0, T). All the above derivatives are continuous in K and 
the limits are uniform in S, y, t on any compact subset of (0, oo) xRx [0, T). 



3 The optimal arbitrage portfolio 

We study the arbitrage from the perspective of the trader, who knows that 
the market is using a misspecified stochastic volatility model to price the 
options. We assume full observation: at every date t, the trader possesses 
the information given by the a-field Tt and knows the deterministic functions 
5", p, a and b. In Section [5] we test the robustness of our results with respect 
to this assumption. 
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To benefit from the market misspecification, our informed trader sets up 
a dynamic self-financing delta and vega-neutral portfolio X t with zero initial 
value, containing, at each date t, a stripe of European call or put options 
with a common expiry date T. In addition, the portfolio contains a quantity 
— S t of stock and some amount B t of cash. 

To denote the quantity of options of each strike, we introduce a pre- 
dictable process (ut)t>o taking values in the space M of signed measures on 
[0,oo) equipped with the total variation norm || • ||y. We refer to [7] for 
technical details of this construction and rigorous definitions of the stochas- 
tic integral and the self-financing portfolio in this setting. Loosely speaking, 
u t (dK) is the quantity of options with strikes between K and K + dK, where 
a negative measure corresponds to a short position. We shall see later that 
for the optimal strategies the measure u is concentrated on a finite number 
of strikes. 

The quantity of options of each strike is continuously readjusted meaning 
that old options are sold to buy new ones with different strikes. In practice, 
this readjustment will of course happen at discrete dates due to transac- 
tion costs and we analyze the performance of our strategies under discrete 
rebalancing in section [5j The portfolio is held until date T* < T. 

Finally, unless explicitly mentioned otherwise, we assume that at all 
times, the total quantity of options of all strikes in the portfolio is equal 
to 1: 



for all t. This position constraint ensures that the profit maximization prob- 
lem is well posed in the presence of arbitrage opportunities (otherwise the 
profit could be increased indefinitely by increasing the size of positions). If 
the portfolio is rebalanced completely at discrete equally spaced dates, and 
the transaction cost per one unit of option is constant, this assumption im- 
plies that independently of the composition of the portfolio, the transaction 
cost paid at each readjustment date is constant, and therefore does not in- 
fluence the optimization problem. This constraint is also natural from the 
point of view of an option exchange or a market maker who wants to struc- 
ture standardized option spreads to satisfy a large number of retail clients, 
because the number of options in such a spread must be fixed. 

Remark 3.1. For an individual trader, who is trying to optimize her option 
portfolio, it may also be natural to impose a margin constraint, that is, a 
constraint on the capital which the trader uses to meet the margin require- 
ments of the exchange. We discuss option trading under margin constraint 
in Section I3TT1 




(7) 



7 



The value of the resulting portfolio is, 

X t = [ P K (S U Y t , t)u t (dK) - 6tS t + B t 



where the subscript K denotes the strike of the option. Together with Lemma 
12.11 our assumptions ensure (see [7J section 4]) that the dynamics of this 
portfolio are, 



, K ^ ,. , , ,^dP K 



dX t =[ I u t (dK)CP K )dt+[ I u t (dK)— ) dS t 



+ u t (dK)^^J dY t - S t d3 t 



where 



d} S'tfcFf Jt&f d 2 f 



dt 2 8S 2 2 dy 2 r dSdy 

To make the portfolio instantaneosly risk-free, we choose, 

C dP K f dP K 

J Ut{dK) -dy- = °' J Ut{dK) -dS~ = 6t 

to eliminate the dY t and dS t terms. We emphasize that these hedge ratios are 
always computed using the market parameters i.e. a, b and p. The portfolio 
dynamics then become 

dX t = J ou t (dK)CP K dt } (8) 

and substituting the equation (jSJ) into this formula and using the vega- 
neutrality, we obtain the risk-free profit from model misspecification: 

dX t = j u t (dK)(C- C)P K dt. (9) 

or, at the liquidation date T*, 

X T * = j j u t {dK)(C - C)P K dt, (10) 

where, 



( r _ r )P K = $M ~ * 2 (Y t )) d 2 P K (bg - b 2 ) d 2 P K 
{ } 2 dS 2 2 8y 2 

+ S t (a t b tPt - a(Y t )b t p t )^ (11) 
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To maximize the arbitrage profit at each date, we therefore need to solve the 
following optimization problem: 



Maximize V t = J u t {dK){C - C)P K (12) 

r r dP K 

subject to / \u t (dK) \ = 1 and / u t {dK)— — = 0. (13) 



dy 

Remark 3.2. It should be pointed out that in this paper the arbitrage portfo- 
lio is required to be instantaneously risk-free and the performance is measured 
in terms of instantaneous arbitrage profit. This corresponds to the standard 
practice of option trading, where the trader is usually only allowed to hold 
delta and vega netral portfolios and the performance is evaluated over short 
time scales. This rules out strategies which are locally risky but have a.s. 
positive terminal pay-off, such as the strategies based on strict local martin- 
gales, which may be admissible under the standard definition of arbitrage. 
Optimality of arbitrage strategies for an equity market under the standard 
definition of arbitrage has recently been studied in |14j . 

The following result shows that a spread of only two options (with strikes 
and weights continuously readjusted) is sufficient to solve the problem (fT2|) . 



Proposition 3.3. The instantaneous arbitrage profit (|12]) is maximized by 

u t (dK) = w]5 K i{dK) - w 2 5 K? (dK), (14) 

where SxidK) denotes the unit point mass at K, (wj, wf) are time- dependent 
optimal weights given by 

dP K i dP K ± 



dP K i , dP K 2 ' * dP K i _|_ dP K 2 

dy dy dy dy 

and (K^,Kf) are time- dependent optimal strikes given by 



^(C - C)P Kl -^(C- C)P K2 
(Kl Kl) = arg max > f*} (15) 



' K 1 ^ 2 &P Kl , dP K2 

dy dy 



Proof. Step 1. We first show that the optimization problem ( fT5i) is well- 
posed, that is, the maximum is attained for two distinct strike values. Let 
F(K l ,K 2 ) denote the function to be optimized in ([TBI) . From the property 
(J6]) and Lemma [27T1 it follows that F is continuous in K\ and K 2 - Let us show 



9 



that for every e > we can find an interval [a, b] such that \F(K X , if 2 )| < e 
for all (if 1 , if 2 ) ^ [a, fe] 2 . We introduce the shorthand notation 

f)pK 

g (K) = —, f(K) = (C-C)P K , 
F{K \ K *) = 9{K 2 )f(K^-g(K^f(K 2 ) 



g(Ki)+g(K* 



Fix e > 0. By Lemma l27Lj there exists N > with |/(if)| < £ for all 

~-\logK\<N <. 



K : K [e , e ]. Let 5 = inf |; gif|<iv f?(if ). Property ([6]) entails that 



<5 > 0. Then we can find M > such that g(if) < eS for all 

if : if ^ [e~ M , e M ] (if the denominator is zero, any M > can be used). It 



gjK T )+gW) 



< 6. In 



follows that for all (if 1 , if 2 ) £ [e^e^] x [ e - M ,e M ] 
the same manner, we can find a rectangle for the second term in F, and, 
taking the square containing both rectangles, we get the desired result. 



Suppose now that for some if 1 and if 2 , F(K l , if 2 ) > 0. Then, by taking e 
sufficiently small, the above argument allows to find a compact set containing 
the maximum of F, and hence the maximum is attained for two strikes which 
cannot be equal since F(K,K) = 0. If, on the other hand, i^if^if 2 ) < 
for all ifi, if 2, then this means that F(K 1 , if 2 ) = and any two strikes can 
be chosen as maximizers of (fT5l). 



Step 2. We now show that the two-point solution suggested by this propo- 
sition is indeed the optimal one. Let uj t be any measure on IR + satisfying the 
constraints ( 1131) . Let Ut = — u>^ be the Jordan decomposition of ut and 

define vf := , , and v7 := , ; . These measures are well-defined due 



to conditions (Q and f lT3|) . From the same conditions, 



The instantaneous arbitrage profit (TT2j) can then be written as 
V t = J{C- C)P K duj+ - J{£ - C)P K du; 

_ j 9 -£dvr j(c - l)P K dvi - / °£dut f(c- t)P K dv^ 
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Then, 



I gdv t I fdvt ~ I gdvf I fdv~ t 



I gdvt + / gdu t 
Jut(dK) {g(K) + jgdv+) sup^ 



g(K 1 )ffdv+-f(K 1 )fgd^ 
g(K 1 )+fgdu+ 



,+ 



< 



I gdvt + / gdv t 
g( K i) [ fdvt - f( K i) I gdvt 

g(K 1 ) + fgdv+ 
Jv+(dK) {gi^) +g(K)} snp K2 



= sup 

Ki 



< sup 



g(Ki)+g(K 2 ) 



= sup 

K X ,K 2 



g{K x ) + fgdv] 
gjK^fjKz) - f(Ki)g(K 2 ) 
giKj+giKi) 



,+ 



Since for the solution (I14p the above sup is attained, this solution is indeed 



The optimal strategy of Proposition 13.31 is adapted, but not necessarily 
predictable. In particular, if the maximizer in (1151) is not unique, the optimal 
strikes may jump wildly between different maximizers leading to strategies 
which are impossible to implement. However, in practice the maximizer is 
usually unique (see in particular the next section), and the following result 
shows that this guarantees the weak continuity, and hence predictability, of 
the optimal strategy. 

Proposition 3.4. Let t* < T , and assume that the coefficients o~ t , b t and p t 
are a.s. continuous at t* and that the maximizer in the optimization problem 
(|l5p at t* is a.s. unique. Then the optimal strikes (K^,Kf) are a.s. contin- 
uous at t* , and the optimal strategy (u t ) is a.s. weakly continuous at t* . 

Proof. To emphasize the dependence on t, we denote by F t (Ki, K%) the func- 
tion being maximized in fjl5[) . By the continuity of o~ t , b t pt, S t and Y t , 
there exists a neighborhood (ti,t 2 ) of t* (which may depend on w 6 Q), in 
which these processes are bounded. Since the convergence in Lemma 12.11 
is uniform on compacts, similarly to the proof of Proposition 13. 'S\ for every 
e > 0, we can find an interval [a, b] (which may depend on u G Q) such that 
\F t {K l ,K 2 )\ < e V(K\K 2 ) £ [a,b] 2 and Vt G (t u t 2 ). This proves that the 
optimal strikes (K^, Kf ) are a.s. bounded in the neighborhood of t*. 

Let (t n ) n >i be a sequence of times converging to t* and {K} , K 2 f) n >i 
the corresponding sequence of optimal strikes. By the above argument, this 
sequence is bounded and a converging subsequence (K\ n ,Kf n )k>i can be 



optimal. 



□ 
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extracted. Let (K l ,K 2 ) be the limit of this converging subsequence. Since 
(Kt n , Kf n ) is a maximizer of F tn ^ and by continuity of F, for arbitrary 
strikes (A; 1 , k 2 ), 

F t ,(K\K 2 ) = lim F Uk (Kl K\ [ ) > lim F t (k 1 , k 2 ) = F t *(k\k 2 ). 

This means that (K 1 , K 2 ) is a maximizer of F t *(-, •), and, by the uniqueness 
assumption, K 1 = K\* and K 2 = Kf*. By the subsequence criterion we then 
deduce that K\ — > K\* and K 2 — > i^ 4 2 * a.s. as t — )■ £*. The weak convergence 
of a;* follows directly. □ 



3.1 Trading under margin constraint 

Solving the optimization problem (fl2j) under the constraint ([7]) amounts to 
finding the optimal (vega-weighted) spread, that is a vega-neutral combina- 
tion of two options, in which the total quantity of options equals one. This 
approach guarantees the well-posedness of the optimization problem and pro- 
vides a new interpretation for commonly traded option contracts such as 
butterflies and risk reversals (see next section). 

From the point of view of an individual option trader, an interesting alter- 
native to ([7]) is the constraint given by the margin requirements of the stock 
exchange. In this section, we shall consider in detail the CBOE minimum 
margins for customer accounts, that is, the margin requirements for retail in- 
vestors imposed by the Chicago Board of Options Exchange 127) . Under 
this set of rules, margin requirements are, by default, calculated individually 
for each option position: 

• Long positions in call or put options must be paid for in full. 

• For a naked short position in a call option, the margin requirement at 
time t is calculated using the following formula |27| . 

M t = P t K + max(aS t - (K - S t )l K>St , f3S t ) := Ff + Xf , 

where P K is the option price, a = 0.15 and (3 = 0.1. 

• The margin requirement for a "long call plus short underlying posi- 
tion" is equal to the call price plus short sale proceeds plus 50% of the 
underlying value. 

• The margin requirement for a "short call plus long underlying" is 50% 
of the underlying value. 
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In addition, some risk offsets are allowed for certain spread options, but we 
shall not take them into account to simplify the treatment. For the same 
reason, we shall assume that all options held by the trader are call options. 
Let u t = oof — uo^ be the Jordan decomposition of oj t . The value of a delta- 
hedged option portfolio can be written as 

X t = B t + j ^P K - S t ^P>uo t {dK) 

= B ' + l{ PK ( l - d ^) +d ^ PK - S <^) 

+ /{- pA '( 1 -^) + ^ (S •- FA ' , }^ (<^A • , 
The margin requirement for this position is 

M t = J \p K + d -^H + l)S^ujt(dK) 

P+(K)u+(dK) + J fr(K)u;(dK) 

with 7 = 0.5. Supposing that the trader disposes of a fixed margin account 
(normalized to 1), the optimization problem (fT2l)- f|T3l becomes 

Maximize V t = J tu t (dK) (C - C) P K 

subject to J 0+(K)u+{dK) + J ^{K)u^{dK) = 1 

r dP K 

and J oj t (dK)-—- = 0. 

This maximization problem can be treated similarly to ( fT2l) - (fT3l) . Assume 
that at time t, the problem 

- C)P Kl - ^{C - C)P K2 
arg sup — -j -5 . (16) 

admits a finite maximizer (K^, Kf). Then the instantaneous arbitrage profit 
is maximized by the two-strike portfolio (1 14)) with 
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However, since (3 + (K) — y as K — > +00, the problem (fT6l) may not admit a 
finite maximizer for some models. For instance, in the Black-Scholes model, 
a simple asymptotic analysis shows that f3 + (K) converges to zero faster than 
as K — » +00, while converges to zero faster than a d g dy and 9 ^ 2 . 
This means that by keeping K 2 fixed and choosing K 1 large enough, the 
instantaneous arbitrage profit ( TT6|) can be made arbitrarily large. In financial 
terms this means that since the margin constraint does not limit the size of 
the long positions, it may be optimal to buy a large number of far out of the 
money options. So the margin constraint alone is not sufficient to guarantee 
the well-posedness of the problem and a position constraint similar to ([7]) is 
necessary as well (in particular, the position constraint limits the transaction 
costs). If both the margin constraint and the position constraint are imposed, 
the optimal portfolio will in general be a combination of three options with 
different strikes. 

While our analysis focused on CBOE margining rules for retail investors, 
a number of other exchanges, for example, CME for large institutional ac- 
counts, use an alternative margining system called Standard Portfolio Analy- 
sis of Risk (SPAN). Under this system, the margin requirement is determined 
as the worst-case loss of the portfolio over a set of 14 market scenarios. In 
this case, position constraints are also necessary to ensure the well-posedness 
of the optimization problem. See |2j for a numerical study of optimal option 
portfolios under SPAN constraints in the binomial model. 



4 The Black-Scholes case 



In this section, we consider the case when the misspecified model is the Black- 
Scholes model, that is, a = b = 0. This means that at each date t the market 
participants price options using the Black-Scholes model with the volatility 
value given by the current instantaneous implied volatility &(Y t ). The process 
Y t still has a stochastic dynamics under the real-world probability P. This 
is different from the setting of [T3], where one also has a = b = 0, that is, 
the instantaneous implied volatility is deterministic. In this section, since we 
do not need to solve partial differential equations, to simplify notation, we 
set at = Y t , that is, the mapping a(-) is the identity mapping. Recall the 
formulas for the derivatives of the call/put option price in the Black-Scholes 
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model (r = 0): 



^ = Sn(di)y/r = Kn{d 2 )\/T. 



da v i/v ™ v ' dS 2 Sa^ 1 

d 2 P n(di)d 2 <9 2 P _ Sn(d 1 )d 1 d 2 ^/r 

dadS a ' <9<r 2 <r 

where cfi ; 2 = ± ^rp, r = T — t, m = log^/i^) and n is the standard 
normal density. 

We first specialize the general optimal solution ( TT4"1) to the Black-Scholes 
case. All the quantities a t ,a t ,b t , pt, S t are, of course, time-dependent, but 
since the optimization is done at a given date, to simplify notation, we omit 
the subscript t. 

Proposition 4.1. The optimal option portfolio maximizing the instanta- 
neous arbitrage profit ( 1121) is described as follows: 



The portfolio consists of a long position in an option with log-moneyness 
mi — Zxa^fr — ^Y- and a short position in an option with log-moneyness 
m 2 = z 2 a^fr — where (zi, z 2 ) is a maximizer of the function 

(>i -z 2 )(zi + z 2 -w) 



f(zi,z 2 ) 



e *?/2 _|_ e z 2 2 /2 



withw=^p^. 

The weights of the two options are chosen to make the portfolio vega- 
neutral. 



We define by P opt the instantaneous arbitrage profit realized by the optimal 
portfolio. 

Proof. Substituting the Black-Scholes values for the derivatives of option 
prices, and making the change of variable z = + ^-rp, the maximization 
problem to solve becomes: 

■naxES \ b -f(4 - 4) - % - *) - ft* (17) 

21,22 n(zi) + n(z 2 j I 2a 2 a 1 



from which the proposition follows directly. □ 

Remark 4.2. A numerical study of the function / in the above proposition 
shows that it admits a unique maximizer in all cases except when w = 0. 
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Remark 4.3. From the form of the functional Q17|) it is clear that one should 
choose options with the largest time to expiry r for which liquid options are 
available. 

Remark 4.4. The gamma term (second derivative with respect to the stock 
price) does not appear in the expression ffTTl) . because the optimization is 
done under the constraint of vega-neutrality, and in the Black-Scholes model 
a portfolio of options with the same maturity is vega-neutral if and only if 
it is gamma neutral (see the formulas for the greeks earlier in this section). 
Therefore in this setting, to arbitrage the misspecification of the volatility 
itself, as opposed to skew or convexity, one needs a portfolio of options with 
more than one maturity. 

The variable z = log j^/fi") _|_ £VpL introduced in the proof is directly related 
to the Delta of the option: for a call option A = N(z). This convenient 
parameterization corresponds to the current practice in Foreign Exchange 
markets of expressing the strike of an option via its delta. Given a weight- 
ing measure Ut(dK) as introduced in section [2} we denote by Ut{dz) the 
corresponding measure in the z-space. We say that a portfolio of options 
is A-symmetric (resp., A-antisymmetric) if u t is symmetric (resp., antisym- 
metric). 

The following result clarifies the role of butterflies and risk reversals for 
volatility arbitrage. 

Proposition 4.5. LetV opt be the instantaneous arbitrage profit (TT2j) realized 



by the optimal strategy of Proposition 



1. Consider a portfolio (RR) described as follows: 

• Ifbr/2 + po/d > 

— buy | units of options with log-moneyness m\ = —Oy/r — 
or, equivalently, delta value iV(— 1) ~ 0.16 

— sell ~ units of options with log-moneyness vri2 = <J\fr — 
or, equivalently, delta value N(l) ~ 0.84. 

• if br/2 + pa jo < buy the portfolio with weights of the opposite 
sign. 

Then the portfolio (RR) is the solution of the maximization problem 
(TL2]) under the additional constraint that it is A-antisymmetric. 

2. Consider a portfolio (BF) described as follows 
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• buy x units of options with log-moneyness m\ = z^o^fr — a 2 T, 
or, equivalently, delta value N(z Q ) w 0.055, where z « 1.6 is a 
universal constant, solution to 

— e 2 = e 2 +1 

• buy Xq units of options with log-moneyness m-i = —z^a^fr — a 2 T, 
or, equivalently, delta value N(—z ) w 0.945 

- 2 

• sell 1 — 2xq units of options with log-moneyness m 3 = — ^ or, 
equivalently, delta value N(0) = \, where the quantity x$ is chosen 
to make the portfolio vega-neutral, that ^tt- « 0.39. 

2(l+e~V 2 ) 

Then, the portfolio (BF) is the solution of the maximization problem 
(I12p under the additional constraint that it is IS.- symmetric. 

3. Define by r p RR the instantaneous arbitrage profit realized by the portfolio 
of part 1 and by V BF that of part 2. Let 

\abr + 2crp\ 

(X = 

\abr + 2a p\ + 2bK oy /r 

where Kq is a universal constant, defined below in the proof, and ap- 
proximately equal to 0.459. Then 

V RR > aV opt and V BF > (1 - a)V opt . 

Remark 4.6. In case of a risk reversal, the long/short decision depends on 
the sign of p + y 2 " rather than on the sign of the correlation itself. This 
means that when the implied volatility surface is flat and p < 0, if b is big 
enough, contrary to conventional wisdom, the strategy of selling downside 
strikes and buying upside strikes will yield a positive profit, whereas the 
opposite strategy will incur a loss. This is because the risk reversal (RR) has 
a positive vomma (second derivative with respect to Black-Scholes volatility). 
Let P x be the price of the risk reversal contract where one buys a call with 
delta value x < \ and sells a call with delta value 1 — x. Then 

QQ. px 

= Sn(d)ard, d = N(l-x)> 0. 

do 1 

Therefore, from equation (ITT]) it is clear that if b is big enough, the effect 
of the misspecification of b will outweigh that of the correlation leading to a 
positive instantaneous profit. 
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Proof. Substituting the Black-Scholes values for the derivatives of option 
prices, and making the usual change of variable, the maximization problem 
(fT2"j) becomes: 

max — — — / z 2 n(z)u t (dz) — Sb(br/2 + pa/a) / zn(z)u) t (dz) (18) 



2a 

subject to / n(z)u) t {dz) = 0, / \u t (dz)\ 



1. In the antisymmetric case it is sufficient to look for a measure on 
(0, oo) solution to 



max —Sb(br/2 + pa/a) / zn(z)u^(dz) 

J(0,oo) 
f 1 

subject to / |w t (<iz)| = -. 

ifO.oo) 2 



'(0,oo) 

The solution to this problem is given by a point mass at the point 
z = argmax zn(z) = 1 with weight — | if br/2 + pa /a > and | if 
br/2 + pa/a < 0. Adding the antisymmetric component of uf, the 
proof of part 1 is complete. 

2. In the symmetric case it is sufficient to look for a measure uf on [0, oo) 
solution to 

Sb 2 



max ■ 



z 2 n(z)uj^(dz) 

[0,oo) 

1 



2a 

subject to / n(z)u^(dz) = 0, / \u^(dz)\ 



[0,oo ) J[0,oo) 



2 



By the same argument as in the proof of Proposition 13. 3[ we get that 
the solution to this problem is given by two point masses at the points 
Z\ and Z2 given by 

n(zi)zln{z 2 ) - n(z 2 )zln(z 1 ) z\ - z\ 

{zi,z 2 ) = argmax -— — = argmax— 3 - T , 

n{z 1 )+n{z 2 ) e^+e-r 

from which we immediately see that Z\ = and z 2 coincides with the 
constant Zq introduced in the statement of the proposition. Adding the 
symmetric part of the measure uf, the proof of part 2 is complete. 

3. To show the last part, observe that the contract (BF) maximizes the 
first term in ( Tl8|) while the contract (RR) maximizes the second term. 
The values of (Tl8]) for the contract (BF) and (RR) are given by 

Sb2 V^^ z ?j2 t>rr _ Sb\br/2 + po-/v\_ e _i 



v bf = ^ v ; c -z?j2 v i 



ctV2tt V27T 
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and therefore 



V RR \abr + 2ap\ 



v bf + v rr \abr + 2ap\+2bK ^ 



with K = e 2 2 . 



Since the maximum of a sum is always no greater than the sum of 
maxima, V opt < V BF + V RR and the proof is complete. 



□ 

Remark 4.7. Let us sum up our findings about the role of butterflies and risk 
reversals for arbitrage: 



• Risk reversals are not optimal and butterflies are not optimal unless 
P = (because a = for this value of p). 

• Nevertheless there exists a universal risk reversal (16-delta risk reversal 
in the language of foreign exchange markets, denoted by RR) and a 
universal butterfly (5.5-delta vega weighted buttefly, denoted by BF) 
such that one can secure at least half of the optimal profit by investing 
into one of these two contracts. Moreover, for each of these contracts, 
a precise estimate of the deviation from optimality is available. 

• Contrary to the optimal portfolio which depends on model parameters 
via the constant w introduced in Proposition I4.1[ the contracts (BF) 
and (RR) are model independent and only the choice of the contract 
(BF or RR) and the sign of the risk reversal may depend on model 
parameters. This means that these contracts are to some extent robust 
to the arbitrageur using a wrong parameter, which is important since 
in practice the volatility of volatility and the correlation are difficult to 
estimate. 

• The instantaneous profit V BF of the contract (BF) does not depend 
on the true value of the correlation p. This means that (BF) is a pure 
convexity trade, and it is clear from the above that it has the highest 
instantaneous profit among all such trades. Therefore, the contract 
(BF) should be used if the sign of the correlation is unknown. 

• When b — > 0, a — > 1, and in this case (RR) is the optimal strategy. 



5 Simulation study in the SABR model 

In this section, we consider the case when, under the pricing probability, the 
underlying asset price follows the stochastic volatility model known as SABR 
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model. This model captures both the volatility of volatility and the corre- 
lation effects and is analytically tractable. The dynamics of the underlying 
asset under Q is 

dS t = a t Sl ( v 7 ! - P 2 dWl + pdW 2 ) (19) 
da t = ~ba t dW 2 (20) 

To further simplify the treatment, we take /3 = 1, and in order to guarantee 
that S is a martingale under the pricing probability [28], we assume that the 
correlation coefficient satisfies p < 0. The true dynamics of the instantaneous 
implied volatility are 

da t = ba t dW 2 , (21) 
and the dynamics of the underlying under the real-world measure are 

dS t = a t 3 t (y/l - pHWl + pdW 2 ). (22) 

The SABR model does not satisfy some of the assumptions of section [2j 
in particular because the volatility is not bounded from below by a positive 
constant. Nevertheless, it provides a simple and tractable framework to 
illustrate the performance of our strategies. A nice feature of the SABR 
model is that rather precise approximate pricing formulas for vanilla options 
are available, see |16| and [23], which can be used to compute the optimal 
strikes of Proposition 13.31 Alternatively, one can directly compute the first 
order correction (with respect to the volatility of volatility parameter b) to 
the Black-Scholes optimal values using perturbation analysis. We now briefly 
outline the main ideas behind this approach. 



First order correction to option price In the SABR model, the call/put 
option price P satisifies the following pricing equation, 

dP S 2 a 2 d 2 P b 2 a 2 d 2 P n 2 ~ d 2 P n 
^ + —dSl + — d^ + Sabp dSda- = ° 

with the appropriate terminal condition. We assume that b is small and look 
for approximate solutions of the form P = P + bP\ + 0(b 2 ). The zero-order 
term P corresponds to the Black Scholes solution: 

dP S 2 a 2 d 2 P 
dt 2 OS 2 



The first order term satisfies the following equation: 

OP, S 2 a 2 d 2 P 1 2 „d 2 P 
dt 2 dS 2 P dSda ' 
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and can be computed explicitly via 



a 2 p(T - 1) d 2 P 
Fl ~ 2 dSd^ 



First order correction to optimal strikes Define 

f = (C-C)P 

where we recall that 



lr ? . n S 2 (a 2 -a 2 )d 2 P (b 2 -b 2 )a 2 d 2 P a/ „ u „ 2 ~. dP 2 
(C - C)P = —— + ^ ; 9 + S{aabp - o 2 bp)- 



dS 2 2 da 2 v r rj dSda 

For the general case, we are not assuming anything on the smallness of the 
parameters of the true model. Now to the first order in b, we can expand / 

as, 

fnfo + bfi 



where, 



and, 



S*{g*-a*) &P b 2~2 d 2p o Qp2 

JO = „ + + Saabp- 



8S 2 2 da 2 r dSda 



S 2 {a 2_^ ) d 2 Pi h 2~ a 2 d 2 Pi #2p2 9 2p2 

h = 2 -dS^ + + S ° ahp dSd-a - Sa P dSd~a 

Define the expression to be maximized from Proposition 13.31 

F _ w 2 f l_ w l f 2 
W 1 + W 2 

where w l = These vegas are also expanded in the same manner: w l 
Wq + bw\. Now to first order in b, we can expand F as, 



F « F := F + b{F 1 + F 3 F ) (23) 



where 



F 
Ft 



and F 3 = 1,0 



Wpfo - Wq/o 
wl+wl 

wl + w 2 

w\ + w\ 
wl + w 2 
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We would like to represent the optimal strike K* which maximizes F as 
the sum of the strike K which maximizes F plus a small correction of order 
of b: 

K* « K + bKi. 
The optimal strike K* is the solution of the equation, 

OF OF ~ d 2 F 

Using the representation (123]) . the fact that K is the maximizer of F and 
keeping only the terms of order up to one in b, we get 

which gives the following correction to the optimal strike: 



§^(K ) + F ^(K ) 

d 2 F a l K \ 



Figure [2] compares this first-order correction to the two optimal strikes 
to the zero-order Black-Scholes value and the precise value computed using 
Hagan's implied volatility asymptotics [TB], for different values of the volatil- 
ity of volatility parameter b. Given that already the difference between the 
zero-order correction and the precise value is not that big, our first order 
correction appears acceptable for reasonably small values of b. 



Numerical results In this section we demonstrate the performance of 
our strategies on a simulated example. We assume that a trader, who is 
aware both about the true and the misspecified parameters, invests in our 
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strategy. We simulate 10000 rung^l of the stock and the volatility during 1 
year under the system specified by (|19)) - (|22p . assuming misspecification of 
the correlation p and the volatility of volatility b but not of the volatility 
itself (a = a). The initial stock value is So = 100 and initial volatility is 
o"o — 1- The interest rate is assumed to be zero. During the simulation, at 
each rebalancing date, the option portfolio is completely liquidated, and a 
new portfolio of options with the same, fixed, time to maturity and desired 
optimal strikes is purchased. 

We illustrate the impact of different characteristics of the strategy on the 
overall performance. In each of the graphs, we plot two groups of curves: 
the cross-marked ones correspond to misspecified model, and the diamond- 
market ones correspond to the setting when the model parameters are equal 
to the true ones. In each group, the three curves correspond (from lower to 
upper) to the 25% quantile, the median and the 75% quantile computed from 
the 10000/100000 simulated trajectories. 

• Impact of correlation misspecification. Figure [3] (left) shows the per- 
formance of portfolios using options with 1 month to maturity. The 
true parameters are p = —.8, b = .2. The misspecified or the market 
parameters are p = —.2, b = .2. 

• Impact of volatility of volatility misspecification. Figure [3] (right) shows 
the performance of portfolios using options with 1 month to maturity. 
The true parameters are p = —.8,b = .2. The misspecified or the 
market parameters are p = —.8,b= .1. 

• Impact of the options' maturity. Figure H] shows the difference in per- 
formance of portfolios using options with 1 month and 6 months to 
maturity. The true parameters are p = —.8,b = .2. The misspeci- 
fied or the market parameters are p = —.2,b= .1. As predicted by our 
analysis in the Black-Scholes case, the performance is better for options 
with longer time to maturity and higher intrinsic value. 

• Impact of the rebalancing frequency. Figure [5] illustrates the effect 
of the total number of rebalancings on portfolio performance. The 
true parameters are p = —.8,6 = .2. The misspecified or the market 
parameters are p = —.2, b = .1. The average profit is roughly the same 
for both the rebalancing frequencies but the variance is higher for lower 
rebalancing frequencies. 

• Impact of trader using different parameters from the true ones. Figure [6] 
shows the performance of portfolios using options with 1 month to ma- 

2 For the figure which includes transaction costs, 100000 simulation runs are done. 
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Figure 3: The evolution of portfolios using options with 1 month to maturity 
for only correlation misspecification (left) and only volatility of volatility 
misspecification (right). 

turity. The misspecified or the market parameters are p = — .2, b = .1. 
The true historical parameters are p = —.8, b = .2 but the trader esti- 
mates these parameters with an error and uses the values p = — .6, b = 
.15 instead. The average profit increases with better estimation of mis- 
specification in parameters. 

• Impact of transaction costs. In all our analysis till now, we have ne- 
glected the transaction costs. In our last test, we include a bid ask-fork 
of 0.45% in implied volatility terms for every option transaction. More 
precisely, if we denote by p mid the option price given by the model 
and by a mid the corresponding implied volatility, the price that our 
trader must pay to buy an option corresponds to the implied volatility 
of a mid + 0.225% and the price she receives for selling an option cor- 
responds to the implied volatility of a mid — 0.225%. This corresponds 
roughly to the bid-ask interval observed for the most liquid options on 
the S&P 500 index. The true parameters are p = —.9,b = .9. The 
misspecified or the market parameters are p = —.1,6 = .1. The evo- 
lution of the portfolio performance with 32 rebalancing dates per year 
using options with 1 month remaining to maturity is plotted in Figure 
[3 With the parameter values that we have chosen, it seems that even 
in the presence of moderate transaction costs, our strategy produces 
a positive pay-off in at least 75% of cases. Whether this will be the 
case in real options markets is a more difficult question, which requires 
an extensive empirical study. We plan to address this issue in future 
research. 
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Figure 4: Impact of the options' maturity: The evolution of portfolios using 
options with 1 month (left) and 6 months to maturity with 256 rebalancing 
dates. 




Figure 5: Impact of the rebalancing frequency: Evolution of the portfolio 
using options with 1 month to maturity and 128 rebalancing dates. Compare 
with Figure HI left graph. 
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Figure 6: Impact of using a different set of parameters: Evolution of the 
portfolio for a different set of parameters than the true ones for options with 
1 month to maturity and 256 rebalancing dates. 




Figure 7: The effect of transaction costs on the performance of the arbitrage 
strategy. 
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A Proof of Lemma 12.1 



First, let us briefly recall some results on fundamental solutions of parabolic 
PDE [21] • Let 

r ( d d\ d 1 xr> / \ d 2 u ^ du , . 

L T *' di u '= at + 2 2. + 2. M*. ^ + ^ *)■ 

x ' i,jr' = l ■> 1=1 

Assumption A.l. There is a positive constant \i such that 

n 

J2 «,,(■<■■ nu, > m 2 v(x,t) g w 1 x [o,t], e e R n . 

i,3=l 

Assumption A. 2. There exists a G (0, 1) such that the coefficients of L are 
bounded and Holder continuious in x with exponent a and Holder continuous 
in t with exponent ~, uniformly with respect to (x,t) in IR n x [0, T\. 

The fundamental solution of the parabolic second-order equation with 
operator L is the function T(x, t, £, T) which satisfies 

£ ^, t, A, ^ ) r(x, f, £, T) = - 0*(* - T), t<T. 



dx dt _ 

Under [A~T1 and |A~2~| the operator L admits a fundamental solution T with 
\DlDF(x,t,t,T)\ < c(T-t)-^exp ("C^f ) (24) 
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where r and s are integers with 2r + s<2,t<T and c, C are positive. 
Consider now the Cauchy problem, 

Lu(x, t) = f(x, t), (x, t) e R n x [0, T), 
u{x,T) = <p{x), xeR n , 

where / is Holder continuous in its arguments, is continuous and these 
function satisfy reasonable growth constraints at infinity. Then the solution 
to this problem can be written as 

u(x,t) = [ T dr [ dST(x,t,Z,T)f(£,T)+ [ d€T{x,U,T)<l>(Z). 



Let us now turn to the proof of Lemma 12.11 We discuss the results for 
the put options. The result for calls follows directly by put-call parity. As 
described in Section [2J the put option price P(S,y,t) solves the PDE (jSJ). 
Let x = log -jl and p(x, y, t) := P(Ke x , y, t). Then p solves the PDE 

!*+Ap = 0, p(x,y,T)=K{l-e x )+- (25) 



where 



1 9 ( d 2 p dp} dp l,r,d 2 p , d 2 p 

a p =o a \ ^4 - \ + + + P ba - 



2 [ dx 2 dx J dy 2 dy 2 dxdy 

The quantities a, b, a, p correspond to market misspecified values, but in this 
appendix we shall omit the tilda and often drop the explicit dependence on 
t and y in model parameters to simplify notation. 

Therefore, the option price can be written as 

p(x, y,t)= [dz [ dvT(x, y, t, z, v, T)K{1 - e z ) + , (26) 



where T(x,y,t, z,v,T) is the fundamental solution of (125|) . Since the coef- 
ficients of A do not depend on x, T(x,y,t, z,v,T) = T(0, y, t, z — x, v , T). 
Coming back to the original variable S, we have 

P(S,y,t) = J dz7(y,^-log^,T)K(l-e z ) + (27) 

where we write 

l{y,t,z,T) := / dvV{0,y,t,z,v,T). 
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Decay of the gamma 



a 2 P 
as 2 



By direct differentiation of (J27J), we get 



d 2 P K K 

= -^i(v,t, log— ,t) 



(28) 



and it follows from (1241) that 



d 2 P 



dS 2 



< 



CK 



S 2 \T-t\~2 



exp 



clog 2 f 



T-t 



(29) 



This proves the required decay properties, and the continuity of also 
follows from 



Decay of the vega ^ We denote U (S, y, t) 



dP{ Q y v,t) and u(x,y,t) 



9p ^Qy'^ ■ Using the regularity of coefficients and local regularity results for 
solutions of parabolic PDE [20, Corollary 2.4.1], we conclude that the deriva- 
tives d ^ig y , Q x Q y i and exist, and therefore the operator ( 1261) may be dif- 



ferentiated term by term with respect to y, producing 



du 



-GO 



d 2 p 
dx 2 



dp 

dx 



(30) 



where 



Ai =A + 



bcr)'— + bb — + a', 
ox ay 



All the primes denote the derivative w.r.t. y and the terminal condition is 
u(S,y,T) = 0, since the original terminal condition is independent of y. 



The right-hand side of ( 13U|) satisfies 
d 2 p dp 



dx 2 dx 



Kj(y,t, -x,T) 



(31) 



so from 



we get 



d 2 p dp 



dx 2 dx 



< 



CK 



(T-ty. 



exp 



cx 



T-t 



(32) 



Let Ti denote the fundamental solution of the parabolic equation with 
the operator appearing in the left-hand side of ( 130]) . Using the estimates 
of the fundamental solution in [2"T] section 4.13] (in particular, the Holder 
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continuity) and the bound (1321) . we can show that the solution to fl30|) is 
given by 

f T f [ dp d 2 P 1 

u(x,y,t) = j dr j ^dzdvT^x^,^ z,v,r)a(v)a (v) <. — - — > (z,v,r). 

Using the boundedness of a and a', the bound on the fundamental solution 
and (1321) . and integrating out the variable v, we get 



Mx,y,t)\<J t dr Jto—p—p exp - j . 

Explicit evaluation of this integral then yields the bound 

\u(x, y,t)\< CKVT -te~^ (33) 
\U(S,y,t)\ < CKVT^te~ C - } ^ (34) 
from which the desired decay properties follow directly. 



Decay of We denote w(x,y,t) — 

differentiate equation (130]) with respect to x: 



du and W(S,y,t) 



<9t 1 \ cte 3 <9x 3 / 



From (ED and (EI 



d 3 p d 2 p 



dx 3 dx 2 



< 



CK 
T-t 



exp 



cx 



T-t 



dU 
dS 



and 



(35) 



Similarly to the previous step we now get: 

f T r f , CK 

\w[x,y,t)\ < dr dz 

' 1 n ~ it J (T-r){r-t) 



-exp 



cz 



c(x — 



T — r 

and explicit evaluation of this integral yields the bounds 



\w(x,y,t)\ < CKe 



' T-t 



\W(S,y,t)\<Cje- — 



77 



(36) 
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Decay of g We denote v(x, y, t) = g, V{S, y, t) - ,, y 
equation (I3"0"j) with respect to y: 



r)(7 



and differentiate 



dv . . d f <9 2 p <9p 

— 1- *4. 2 W = -crcr — | 



9* 



cte 2 



9p\ // / , V' 
— > — a u — (pba) w 

ox ) 



where 



A 2 =Ai + bb'— + (pba)'— + a' + {bbj. 
ay ox 



Once again, from (|3T|) and ()2 



9 


J 9 2 p 


dp 1 


9y 


\ dx 2 


9x / 



< 



T-t 



cxp 



car 



Using this bound together with (136|) and (133|) and preceding as in previous 
steps, we complete the proof. 
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